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Abstract 

Molecules of a nematic liquid crystal respond to an applied magnetic field 
by reorienting themselves in the direction of the field. Since the dielectric 
anisotropy of a nematic is small, it takes relatively large fields to elicit a 
significant liquid crystal response. The interaction may be enhanced in col- 
loidal suspensions of ferromagnetic particles in a liquid crystalline matrix — 
ferronematics — as proposed by Brochard and de Gennes in 1970. The ability 
of these particles to align with the field and, simultaneously, cause reorien- 
tation of the nematic molecules, greatly increases the magnetic response of 
the mixture. Essentially the particles provide an easy axis of magnetization 
that interacts with the liquid crystal via surface anchoring. 

We derive an expression for the effective energy of ferronematic in the 
dilute limit, that is, when the number of particles tends to infinity while 
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their total volume fraction tends to zero. The total energy of the mixture 
is assumed to be the sum of the bulk elastic liquid crystal contribution, the 
anchoring energy of the liquid crystal on the surfaces of the particles, and 
the magnetic energy of interaction between the particles and the applied 
magnetic field. The homogenized limiting ferronematic energy is obtained 
rigorously using a variational approach. It generalizes formal expressions 
previously reported in a physical literature. 

Keywords: Ferronematics, liquid crystal, homogenization 



1. Introduction 

The study of magnetic particle suspensions in a liquid crystalline matrix 
was initiated with the theoretical article by Brochard and de Gennes 

Q (July, 

1970), and the experimental work carried out by Rault, Cladis and Burger [2|, 
(June, 1970). 0. The underlying mechanism behind a ferronematic system 
is a mechanical coupling between the nematic molecules and the magnetic 
particles, mostly realized by the surface anchoring energy. 

Molecules of nematic liquid crystals have positive magnetic susceptibility, 
so they tend to align themselves in the direction of an applied magnetic field. 
However, since this magnetic susceptibility is small — of order 10~^ — it takes 
large fields, about 10^ Oe, to elicit a significant response. Brochard and 
de Gennes argue that the addition of paramagnetic ions to the system is 
not an efficient way to increase the magnetic susceptibility constant, since it 
would require a concentration of paramagnetic ions above n = 10^° ions per 
cm^. The latter is the limiting value that cannot be exceeded in order to 
guarantee the preservation of the liquid crystal properties of the system. 

The focus of research turned to suspensions of large ferromagnetic parti- 
cles in the nematic matrix. Brochard and de Gennes identified the two key 
properties of such systems: strength of the mechanical coupling and stability 
of the suspension. The former guarantees that the effect of the magnetic 
field on the liquid crystal, acted through the magnetic particles, is ability 
to control the nematic texture. The latter property sets a limit on the size 
and concentration of particles to prevent clustering. The numbers arrived 
at from theoretical considerations set the particle length I > 0.5 x lO^^/xm, 
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underway. 



2 



and a ratio ^ ^ 10, where d denotes the diameter of the particle. The theo- 
retical prediction on particle volume fraction was not to overcome the value 

/ = 10-^. 

In their experiments, Rault, Cladis and Burger |2| chose mono-domain 
particles of 7Fe20|, of 0.35/im long (1) by 0.04 /xm in diameter (d). The sat- 
uration magnetization is 384 gauss with the easy axis parallel to the long axis 
of the grain. Grains of these dimensions satisfy the criterion for mechanical 
coupling to the nematic liquid and also for mechanical rotation, as opposed 
to rotation of magnetization inside the grain, in a reversed field. Typical 
grain concentrations were of the order of 2 x 10^^ grains/cm^, which corre- 
sponds to / ~ 1.4 X 10~^, well within the theoretical prediction by Brochard 
and de Gennes. For this physical parameters, Rault, Cladis and Burger state 
jif: The ferronematic appeared to be very stable in the nematic-isotropic 
phases showing very httle tendency to agglomerate. However, if a high field 
(1 kg) is applied to the sample in the isotropic phase, upon returning it to 
the nematic phase, we have observed long chains of grains about 50fim. Both 
works assert that distortions of the nematic pattern in magnetic suspensions 
occur at very low fields: magnetizations range in the order of 0.1 to 1 gauss, 
instead of the values 10""^ to 10^^ of pure nematic liquid crystals, with a 
typical coupling gain of order 10^. 

Central to the understanding of the nematic-magnetic coupling is the 
question of how the grains align in the nematic. Brochard and de Gennes 
postulated strong anchoring of nematic molecules along the magnetic mo- 
ment, assumed to coincide with the direction of the particle axis. The effect 
of the grain magnetic field results from the anisotropy of the field around the 
grain, present even in the case of a spherical grain, resulting in a preferential 
direction for the magnetic moment in the nematic phase. This effect turns 
out to be small for small grains, with the magnetic moment causing only a 
local disruption of the nematic alignment. 

In their experimental work, Chen and Amer [jij, 1983] used particle coat- 
ing that yields homeotropic anchoring of the liquid crystal on the magnetic 
grain to synthesize stable ferronematic systems. Although the length and 
aspect ratio of the particles, 0.5 /im and 7 : 1, respectively, are compatible 
with those considered in the previous works, the earlier theory assuming rigid 
parallel anchoring was found to be not applicable to the homeotropic case. 
The question of the surface anchoring and its implication on the relative 
orientation of m and n gave rise to an intense experimental and theoretical 
research activity spanning over three decades. In jsj, joj), the authors 
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showed that the rigid anchoring approximation, m||n might be used only 
if the condition ^ >> 1 holds, where W represents the surface energy 
density, and K denotes a typical Frank constant. A calculation for MBBA 
data, with K = 5x 10"^^, and IQ-^ <W < 10'^ and d = 0.07/im gives 
lQ-2 ^ wd, ^ 10^^, showing a finite surface energy of the system. 

Assuming soft liquid crystal surface anchoring, Burylov and Raikher ([3], 
1995) proposed a macroscopic free energy density of the form 

F = -{Ki{dwny + K2{curln- n)^ + ^3(11 x curlnf} - -Xa(n ■ H)^ 

:M.fi^.HH{l^ynf^{^)(n.mr (1) 

Here / represents the volume fraction of the particles, Xa the anisotropic 
part of the diamagnetic susceptibility of nematic, and the positive constants 
u and Ms denote the particle volume and the saturation magnetization, re- 
spectively. In the last term, W represents the strength of the surface energy 
and y4 = 1 — 3 cos^ a characterizes the type of anchoring, with a denoting 
the easy-angle orientation of the nematic on the particle surface. 

The macroscopic free energy ([1]) has been investigated in theoretical and 
experimental works involving orientational transitions in ferronematic states 
[sl, [sj, j^. In particular, ^] presents a nonlinear modification of the Rapini- 
Papoular energy that predicts a first order Fredericks transition. In and 
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Kopcansky et al. use the modified theory to determine threshold fields 
in ferronematic transitions under combined electric and magnetic fields. In 



12l |. the authors report on experimental studies of structural transitions in 
ferronematic subject to electric and magnetic field, with the matrix consisting 
of 8CB and 6CHBT liquid crystals, respectively. While in both cases the 
anchoring was determined as soft, it was found that n ± m in the first case, 
and n||m in the second. So, it was established then that both, parallel and 
perpendicular anchoring may occur depending on the properties of the matrix 
(which, in turn, reflects the properties of the particle coating). Zadorozhnii 
et al. [l3| provide a comprehensive analysis of the director — a unit vector in 
the direction of the preferred molecular alignment — switching for small and 
large values of the applied field in a nematic liquid crystal cell subject to 
homeotropic boundary conditions at the cell and particle walls. They show 
that the threshold field depends on the anchoring strength of the director on 
the particle surface. 
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Note that a closely related set of models ij]- 15| exists for suspensions 
of ferroelectric nanoparticles in a nematic liquid crystalline matrix. The 
mechanical coupling between the particles and the nematic is still governed 
by the surface anchoring, but the particles interact with an electric and not 
a magnetic field. 

In this work, we rigorously derive an expression for the effective ferrone- 
matic energy that reduces to the models described above under appropriate 
limits. We consider a collection of spheroidal particles with, fixed, randomly 
distributed locations in the matrix, and with magnetic moment pointing in 
the direction of an easy axis. The particles are taken as rotations and trans- 
lations of the same spheroidal particle, located at the origin. We model 
the liquid crystalline matrix according to Ericksen's theory of nematics with 
variable degree of orientation. In this theory, the state of a liquid crystal 
is described by a vector u(x) whose direction gives the average molecular 
alignment at the point x, while its magnitude |u(x)| — degree of orienta- 
tion — describes the quality of the alignment. Assuming that the Frank 
elastic constants are equal, the bulk liquid crystal energy has the form of 
the Ginzburg-Landau energy for u. We assume soft anchoring of the liquid 
crystal molecules on the surfaces of ferromagnetic particles as represented by 
the Rapini-Papoular energy term. The surface energy contribution can be 
either positive or negative depending on whether parallel or perpendicular 
alignment of nematic molecules is preferred on particles surfaces. It turns 
out that the case when the surface energy is negative is the most challenging 
to analyze. 

Mathematically, we consider a family of energy functionals, J-'^, parametrized 
by a quantity e > that characterizes the geometry of the system, specifi- 
cally, the size of the particles and the inter-particle distance. The system is 
assumed to be dilute, that is the volume fraction of the particles tends to 
in the limit e — )■ 0. The parameter scalings of the model that give the rela- 
tive contribution of the different components of the energy are formulated in 
terms of e as well. The choice of scalings guarantees that the limiting contri- 
butions of the bulk and surface energies, as well as the energy of interaction 
between the particles and the applied magnetic field are of order 0(1). We 
show that for the same parametric regime the contribution from the energy 
of magnetic interaction between the particles is o(l) in e. This is consistent 
with the experimental observations that characterize dilute small particles 
systems in the absence of clustering. 

We study the variational limit of the family of energies {J^e} as e — )■ 0. 
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The limiting functional {/"o} represents the effective, or homogenized energy 
of the system. Here the convergence is understood in the sense that the 
sequence of minimizers {u^} of {J^;} converges to a minimizer u of {J^o} in 
an appropriate functional space. The effective energy provides a benchmark 
for comparison with the formal expression for ferronematic energy functional 
0] given in (DQ). 

The homogenized energy f lTB]) is more general than ([T]) as it is obtained 
under less restrictive assumptions. The interaction between the liquid crystal 
and the particles is due to surface anchoring and is represented by the matrix 
A in (fTSjl that encodes the information on the shape and size of the particles, 
their locations, and their orientation with respect to a fixed frame. Likewise, 
the effective magnetic moment M in (1T5]) that couples the particles to the 
external magnetic field depends on the spatial and orientational distributions 
of the particles. For the high- aspect-ratio, needle-like particles the coupling 
terms reduce to their counterparts in ([1]). 

2. Background 

Given the domain Q C let G flhe an arbitrary collection of subsets 
of fl such that Pi fl Pj = for every i j where i,j = 1, . . . , n. Suppose 
that the region Q\ Uj Pi is occupied by a nematic liquid crystal and that for 
each i = 1, . . . , n the region Pi corresponds to a hard ferromagnetic particle 
embedded in the nematic matrix. 

We will consider the liquid crystal configurations that can be described 
by the Ericksen's theory for nematics with variable degree of orientation; we 
will neglect all flow effects and assume that all elastic constants are equal. 
Further, we will use the phenomenological Rapini-Papoular term in order 
to approximate the liquid crystal/ferromagnetic surface energy. Then the 
elastic energy of the liquid crystal is given by 

J^ic-= [ {K\Vu\^ + Wi\u\)) dV + q [ {u,uYda, 

where i^' > is the elastic constant, g G R is the strength of the surface term, 
W is the bulk free energy of the undistorted state, and u is the outward unit 
normal vector to dPi. 

Suppose that ferromagnetic particles are sufficiently small so that for 
every i = 1 , . . . N an i—th particle can be characterized by a magnetization 
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vector rrij pointing in the direction of an easy axis of the particle. In order 
to derive the expression for the magnetostatic contribution /™ to the free 
energy density of what is effectively a diamagnetic matrix interspersed with 



the ferromagnetic particles, we follow 16|. We have that 

'df 



dH 



= -B , (2) 

T 



where H and B are the magnetic field and the magnetic induction, respec- 
tively (cf. eq. (39.1) in (Toj) and the derivative is taken holding the tempera- 
ture T fixed. Assuming that M denotes the magnetic moment, the induction 
is given by 

B = /io(H + M), (3) 

where /iq is the magnetic permeability of vacuum. 

Suppose that the magnetic moment of the material can be written as 

M = m + xH , (4) 

and the material can exhibit both the spontaneous magnetization m (an 
independent thermodynamic variable) and the magnetization induced by the 
field (we assume it to be proportional to the field). The tensor x is the 
magnetic susceptibility; it is generally small in diamagnetics, but it can be 
large in soft ferromagnetic bodies. In what follows, we will set m = in the 
liquid crystal while we will set x = in hard ferromagnetics. 

Substituting ([3]) and (jl]) into ([2]) and integrating with respect to the field, 
we obtain 

r (m, H) = r (m, 0) - /io(m, H) - (5) 

Here /i = /io(I + x) is fhe magnetic permeability tensor. Note that the 
energy /'"(m, 0) accounts for both the exchange and anisotropy energies for 
a ferromagnetic body. We will ignore this splitting since we consider single- 
domain particles. 

The expression ([5]) can be adjusted further by excluding the energy of 
the external field that would otherwise be created by the same sources in 
vacuum. 

Let the fields H and h solve the (different) sets of Maxwell's equations 
under the same boundary conditions at infinity in the presence and in the 
absence of the material, respectively. Then h is the magnetic field in vacuum 
when there is no magnetizing body (cf. eq. (32.1) in |l6l|). 
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Since the free energy of the field h is 

2 



fio[h 

c 

the adjusted free energy can be written as 



R3 ^ 



By rearranging terms, using Maxwell's equations, and integrating, one can 
show 161 that 

B 'J 

^"^^ L (^'" + ^(H,B)-^(M,h)^ dV. (7) 

This equation can be simplified by taking into account to obtain 

:F^ = -^[ ((m,H) + (m,h)+x(H,h)) dl-, (8) 

where we dropped the tilde for convenience. 

In a hard ferromagnetic material, the magnetic susceptibility x = 0- By 
denoting the demagnetizing field by Hj = H— h the equation (jH]) reduces to 

•^''" = -V / ((m,H,)+2(m,h))dV^, 

— this is the sum of the magnetostatic and the Zeeman energies. Further, 
Hj vanishes as a; — ?■ oo and it satisfies the same set of Maxwell's equations 
as H. 

If the material is diamagnetic, then m = and x is small enough so that 
the magnetic field is essentially unperturbed by the presence of magnetizing 
body. We conclude that 

r^ = -^ I x(h,h)d\/. 
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R3 



which is the standard form of the free energy for the diamagnetic bodies. 

Now we establish the expressions for the magnetic free energy in various 
components of the composite. Suppose for now that the external field h is 
constant. 
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Using the same notation as above, the energy of interaction between the 
magnetic field and the (diamagnetic) hquid crystal (cf. 18|) is given by 



ic ■— ~~ 



2 



X,c(H,h) dV. 



n\UiPi 



The magnetic susceptibihty tensor xic can be approximated as 

Xa|U| / U U 1. 



Xlc 



^exp 



U U 

u u 



Here Xa = X\\ ~ X± is the rescaled diamagnetic anisotropy and 

X = {X\\ + 2x±) /3, 

is the average susceptibility. The scaling factor Sexp is the value of the uniaxial 
order parameter |u| when the measurements of the susceptibility were taken, 
and it refiects the hysteresis behavior of the magnetic loading experiments. 
We point out that in a nematic x\\ < and < Xa ^ Ixl [ISl- The 
smallness of Xa/ X is the basis for assuming that the effect of the liquid crystal 
on the magnetic field is weak 17 . 

By setting x = in ([8]) , the free energy of the hard ferromagnetic particles 

is 



N „ 

ifY. / {K,H) + (m„h)}rfy. 



By solving the Maxwell's equations of magnetostatics, we find that the total 
field H is given by 

H = -V(/., 

where the magnetic potential satisfies the equations 

A0 = 0, in U Pi, 
div(/i,,V0) = 0, in {UPiY. 

The boundary conditions are 



0, 



dPi 



for every i = 1, . . . , N and 



V0 = h, 
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at infinity. Here the magnetic permeability tensor /i = /i/c = /io (I + Xic) in 
the hquid crystal and /i = /xqI in the ferromagnetic particles. 

The equilibrium configuration of the composite can be found by minimiz- 
ing the functional 

with respect to u and mj. 



3. Formulation of the problem 

Suppose that the positions and orientations of prolate spheroidal particles 
are fixed and distributed randomly in the matrix, the spontaneous magnetic 
moments of the ferromagnetic particles are parallel to their long axes, and 

Xa = 0. 

Consider the family of energy functionals J-'^ 

jr [u] = /" {\Vu\^ + W{\n\)}dV + ge [ {u^vf da 



/ {(m„H,) + 2(m„h,)}dV^, 

JR3 



where e > is a small parameter related to the geometry of the system. Here 

m, = mf, xGPf , . 

0, x G \ upf . ^ > 

and for simplicity, we set W{t) = (1 — t^)^. The magnetic field is given by 
he = Ih^l = constant, = —Vip, (11) 

with 

Av? = 0, X e R3, 

[-/.o|^ + (mf,z.)]|^^^=0, xedV^. ^12) 

We assume that for a prescribed U G C^{Q, R'^), 

u, = U, on dn. (13) 

For each e > 0, we denote by the minimizer of (Q. We study the limiting 
energy and the behavior of minimizers of J-'^ as e — > 0. 
We make the following assumptions: 
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1. The ferromagnetic particles consist of a family of prolate spheroids 

= + e"RlV, i = l,...A^e, where G denotes a particle 
center and V is a. reference spheroid with the long axis parallel to the 
z-coordinate axis, and i?f is a rotation. 

2. Given positive numbers < d < D, the distance between particles 
|xf - x^l G [de,De], for all, < < N^. Thus < Ne~^ for some 

> uniformly in e. 

3. |mf| = rrie = Vol{Vf)me'^^ , = he'^^, and Qs = ge"' , where m, h, and 
g are given constants. 

4. The parameters a,/9i,/92, and 7 satisfy 

1 < a < 2, 6a + 2/3i > 9, 
/32 + /3i = 3 — 6a, 7 = 3 — 2a. 

5. The matrix- valued functions 



(14) 



(x) = e^mVof (P) x^) i?fz 



(15) 



converge in the sense of distributions to A, M G L°°(f2), respectively, 
where A : — )■ M^^^ and M : — )• M^. Here z is a unit vector in 
the direction of 2;-axis. 

Remark. Note that the total volume of the particles satisfies Vol (UPf) = 
O (e:'^'-""^-'), so that the homogenization problem for ^ corresponds to a 
dilute limit when limVol(U'Pf) — > as e — ?■ 0. The scalings on and g^ 
guarantee that the magnetic interaction between the applied field and the 
particles, the Ginzburg-Landau energy, and the surface energy are all 0(1) 
while the magnetic interactions between the particles are of order o(l) and, 
therefore, can be neglected. 

Our principal goal is to prove the following 

Theorem 1. Suppose that the assumptions 1-5 hold. Then the sequence of 
minimizers {u^j^^Q of the functional J^el^e] converges in the sense of (d^) 
to a minimizer of the functional 

J-Q^u] = f [|Vu|^+(l-|u|2)%(Au,u)-2(h,M)] dy, (16) 

where A and M are as defined in assumption 5. 
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The matrix A and the vector M that appear in the statement of Theorem 
[U describe the homogenized hquid crystal/ferromagnetic particles interaction 
and the effective magnetization density, respectively. 



4. Main Results 

We prove Theorem [1] in several steps as outlined below. 

4.I. Liquid Crystal Energy 

First, we consider the energy (|9]) without the magnetic terms, that is 

£e[u]= [ {\Vu\'' + Wi\u\)}dV + ge [ {vi.pfda. (17) 

For each small e > 0, we let be the minimizer of flTTl) subject to the 
Dirichlet boundary condition = U on dVt. 

We want to find the limiting functional of the family as e — )■ 0. Al- 
though our approach is developed for the prolate spheroidal particles, it can 
be easily extended to particles of arbitrary convex shapes. The method is 
based on the procedure developed in i^] for the case of spheres. 

^.i.i. Compactness 

We first observe that the restriction of U to the domain VLg, = VL \ VJVf is 
an admissible function. Indeed, 

£,[\]]= [ {\VV\^ + W{\V\)}dV + g, [ {V,u)da 

< [ {\VV\^ + W{\V\)}dV + ge [ {V,u)da ^^^^ 

<C (1 + geN,\dV\e^'') < C (1 + gN\dV\) < C, 

where C is a generic positive constant. Consequently, 

£e[n,]<£,[V]<C. (19) 

That is, ^^^[ue] is uniformly bounded in e. 

The following lemma is needed towards the proof of compactness of the 
sequence {ug} of energy minimizers of ([Tj 
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Lemma 1. Let V denote a prolate spheroid in with minor and major axes 
A and B, respectively. Let V D V represent the prolate spheroid homothetic 
to V with axes ^ = ^ > 2. Then 



/• , ,2 , 352(1 + A) [ ,2 

JdV ^ Jv\V 



A 



7 A Jvyp 



dV 
\u\UV. 



(20) 



Proof. Suppose that the center of the spheroid V is at the origin and its long 
axis is oriented along 2;-axis. We introduce the coordinates 

X = p sincj) cos 9 , y = psincpsmO, z = A^^ Bp cos cj), 

then the volume element is given by dV = A'^Bp^ siiKp dp d6 d(j) and p = C 
defines a prolate spheroid with axes C and BC/A with the surface area 



element da = A ^C'^ sin (f)\/B'^ siii^ (p + A'^ cos"^ (p d9 dep. We start with the 
relation 



u{t,(p,6)— / Upd p, where t G [A, A]. 



(21) 



Let A > be fixed. Taking the square of (|2T1) and applying Young's inequality 
gives 



\uf{A, 0, 9) = \uf{t, 0, 9) - 2u(t, 



Updp + 



Updp 



<(1 + A) 
Further, by Holder's inequality 



Updp + ( 1 



(22) 



\u\\t,(P,9). 



/ Updp < / \upfp'^dp / p "'dp 

J A J A J A 
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We multiply both sides of the inequality f l22p by the determinant of the 
Jacobian, integrate inV \ V, and use the fact that |up|^ < 3-B^A~^| Vu|^: 



B r 
1 



2n rA 



|u| {A,(j),6)p sm(j)dpd9 1 

Jo J A 



<—^ ^ / / / \up\^ A- 'Bp^ sin (j) dp d9d(p 



1 



JO JA 

TT r2-K i-A (23) 
|2 /1-1d„2, 



+ (l + -jy j J \u\' A- ' Bp^ sin (p dp de, 
B^(a^- A' 



<(1 + A) ^— [ \VufdV+(l + \][^ \ufdV. 

A J.p\^p \ A J Jfyp 

At the same time 

r-TV f2Tt 



A 



|u| {A,(f),9)p sin (j) dp d9 d(j) 
Jo JA 

r '''' ^ ' 

y|3 _ /-vr /•27r 



, (24) 

> -— / \uf{A,^,9)asin^JB'^sin'^^ + A'^cos'^(pd9d(p 

3^ Jo Jo 



A^-A^ 
3A^ 



lup da . 



dV 



Combining ( l23l) with ( 1241) and using the fact that A > 2A, we obtain 

,2, 352(1 + A) f ,o / 1\ 3^2 f , ,2,.. 



dv ^ Jv\v \ ^JA^ — A^ Jv\v 

3^2(1 + A) /■ ,2 / 1 A 24^2 /■ , ,2 

A Jv\v \ Ay 7^3 Jv\r 



□ 



Next, we use the previous lemma to estimate the surface energy contri- 
bution in (fTTIl in terms of the L2— norms of u and Vu. 
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Lemma 2. Let e > 0, X > be as in LemmaUi Then 



g, I {u-uyda < C(l + A) 



+A-^ / \u\^dV 



e I \Vu\'^dV 



(25) 



for any admissible function u, where the constant C is independent of e. 

Proof. Let C denote a generic constant independent of e. Setting A = e"a, 
B = and A = de/2, we apply Lemma [1] to the surface integral term 



(u ■ da < / I u| ^ (icT 

dPf JdVf 



< e 



.36^(1 + A) 



\Vu\'dV + e 



2a-3 



1 + T 



1\ 192a2 



A/ 7rf3 



\n\'dV 



Vf\Pf 



< C(l + A) 

JVf\Pf Jvi\vf 
Then, since Qe = ge^~'^°' and 1 < a < 2, we have 

ge (u ■ z/)^ da = gey2 (u ■ i^f da 
<C(1 + A)^ 



i=l 



e / \VuydV + \-'^ / \uydV 



(26) 



<C(1 + A) 



e / iVupdl^ + A"^ / \MVdV 



□ 



We are now in the position to prove the following theorem 



Theorem 2. If a sequence of admissible functions {ug} satisfies [u^] < M 
for some constant M > uniformly in e, then there exists a constant M > 



such that 1 1 Up 



\m{Q.\yjiV<') ^ ^ uniformly in e. 
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Proof. Suppose that {u^} satisfies [u^] < M for some constant M > 
uniformly in e. Using Lemma [2] with A = 1, the assumption on W{t), and 
Holder's inequality, we have 

[ {\Vuf + \uf} dV < M + 2 [ \ufdV+\g,\ [ {n ■ uf da 



- \n\\JVt\ < Ce 



WuVdV 



n\uvf 



n\uvf 



+ C 
< Ce 

in\uvr 



\urdV + Mi 



WurdV 



\u\Uv] +Mi, 



where Mi > is a constant independent of e. Let e be small enough so that 
Ce < |. Then 



n\uvf 



iVul + lul 



dV < Mo 



1 + 



\ufdV 



uniformly in e for some constant M2 > 0. Using the same arguments as in [2( 
we conclude that there exists a constant M > such that ||ue||j:^i^j^^u.pE'j < 

M uniformly in e. □ 

Remark 1. Note that the proof of Theorem is trivial if g > when the 
boundary term is nonnegative. 

Due to our assumptions on the distributions and the sizes of the spheroids 
the domains in the sequence Q \ UPf are strongly connected [211] that is, 
for every function u G {Q \ UPf ), there exists an extension u G H^[fl) 
such that 

where C > is independent of e. Note that a sufficient condition for f l27|) is 
the existence of a " security layer" around each particle having thickness com- 
parable with the diameter of the particle as e — )■ 22| . It follows that there 
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exists a sequence {u^} of extended minimizers that is uniformly bounded in 
H^{Q) and, up to a subsequence, converges to some Uq weakly in H^[fl) and 
strongly in L'^{Q). Thus 



/ 



\u^-UofdV^O, (28) 



in\urf 

as £ — >■ 0. Further, by trace theorem, 

(uo-U)|an = 0. 

In order to identify the limiting functional and to demonstrate that Uq is 
its minimizer, we now prove 

Theorem 3. Suppose that 

S[u] := J [l Vu|' + (1 - |u|')' + {Au, u)] dV , (29) 

for every u G H^{fl). Given w e C°^{fl), there exists a sequence {w^} C 
H^{fl) such that 

S,[w']^£[w], (30) 

when £ — )■ 0. 

Proof. We begin by constructing a test function. Let w e C°°(0) and set 
^^:=y, + ^^ = ^ + Y^(nt-w)<pO^^) , (31) 

i ^ ' 

where k e (1, a), the function e C°°(M"'") satisfies 

'^^^^ \ 0, if i > 1 , 
For every i — 1, . . . ,N^, the function is a solution of the following problem 



Auf - ^ (uf - w,) = 0, in B,. (x|)\Pf , 

^ + g^{nl,iy)iy = 0, ondVf, (32) 

uf = Wj, when |x| = e'^ , 



where Wj = w(x|) 
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To understand the behavior of a solution to (152]) . for a fixed i G {!,..., A^g}, 
we rescale the lengths by the characteristic size of the particle: y = £:^°(x — 
xf) and set uf (y) := uf (xf + e^y) — Wj. Then 

Auf - uf = 0, in B,.-^{0)\Vi , 

^ + ge''-"{ul + Wi,u)iy = 0, on dVi , (33) 
uf = 0, when |y| = e*^"" , 

where the spheroid Vi = e~°'Vl is centered at the origin. Note that uf is a 
critical point of the functional 

Ei[u]:= [ [\Vuf + \n\'^] dV + ge^-'' [ {n + Wi,ufda, (34) 

where u G Hq {BR{0)\Vi) and R = e'^~°'. We can assume that uf is a global 
minimizer of El over {B^{0)\Vi) once we prove the following 

Lemma 3. The minj|^i(-^^(o)\Pi) ^1 is attained and the minimizer satisfies 

dV < Ce^-^" , (35) 



2 

ulf da < Ce^-^"" . (37) 



\u^f dV < Ce^-^" , (36) 



'dVi 

Proof. 1. Boundedness from above. Since u = is in {Bii{0)\Vi), 

min Ei < El [0] = ge^''' [ (w(xf), da < 1 , (38) 
j?i(i?fl(o)\n) Javi 

when e is sufficiently small. 

2. Boundedness from below. When g > 0, the result is automatic as the 
functional El is nonnegative. Suppose that g < 0. Let e > be small enough 
so that Vi C Bji{0) and choose u G (-B/j(0)) such that the support of u 
is contained in Bfi{0). Following the same line of reasoning as in the proof of 
Lemma [Hand switching to spherical coordinates with 2;— axis along the long 
axis of the spheroid Vi, we have 



u(p(0),^,0) = - r nr{r,e,<P)dr 
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where 

p{4>) = (39) 

[b'^sm (j) + a^cos^^j ^ 

is the equation of the spheroid. By Holder's inequahty 

R \^ 

Wr{r,6 ,(j)) dr\ < I |ur(r, ^, </))|^r^ / r'"^ dr 
pW / Jpi4>) Jp{4>) 

< — — — - / \\ir{r,6,(j))^r'^ dr , 







■f 


Ur (r. 






1 






f 








'pW 



then 

\u{p{4>),e,4>)f <^ I |u,(r,e,0)|Vrfr. (40) 

For the prolate spheroid with long axis in the direction of 2;-axis, the element 
of the surface area is given by 

da = (p^ + pI) V sin (pdOdcf). (41) 
Multiplying fHOl) by the Jacobian and integrating, we obtain 



u(p(0) , 6, 0) I' (p2 + '/'p sin d9 , 



^0 



< 



r r r \ur{r,e,(l))\\p^ + pl^^'^sincl^drde^ 

Jo Jo Jp((l>) 



then 



/ \ufda< max (p^ + pj)^^^ / \Vuf dV 



< C I iVul^dV. 



(42) 



where the constant C > depends only on Pj. 
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Wj, da 



Using fl42l) we obtain the following estimate 

Ei[u]= [ [\Vuf + \uf] dV + ge'^-'' [ (u + 

>[ [\Vuf + \uf] dV - 2\g\e^"' [ {u,ufda 

+ [ {wi,ufda ={1-06^-") [ \VufdV (43) 

+ [ \uf dV - 2\g\e^-'' [ {wi,ufda 

>\f [|Vu|' + |u|^] dV-l, 

^ Jbr{o)\p, 

when e is sufficiently small uniformly in u. It follows that 

El[u]>-1, 

for the same values of e. Since {BR{0)\Vi) is dense in Hq {BR{0)\Pi), 
the inequalities (gS]) and (jH]) hold for all u E {BR{0)\Vi). 

3. Existence of a minimizer. Suppose that {u^} C {Bji{0)\Vi) is a 
minimizing sequence for El. For a sufficiently small e, from (155]) and ( I43p we 
can assume that 

/ [|Vu|V |u|^] dV <2, (44) 

uniformly in k. Then, up to a subsequence, {u^,} converges weakly in the 
space i^Q (i?j?(0)\Pi) to a that minimizes El by the lower semicontinuity 
of (IM]) and the trace theorem. 

4. Properties of the minimizer. In this part of the proof, C denotes 
various constants that depend on Vi and Wj only. Multiplying the equation 

by uf and integrating by parts over Bji{0)\Vi, we have 

/ + \utf] dV 

JBn(Q)\Vi ^^^^ 



-ge^ " / (uf + Wj,z/) (ulu) da. 



1 1 
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From f l43p and Holder's inequality it follows that 

[ {ul + Wi,u){ul,u) da < [ {ulvfda + c(! {ulufda 

JdVi JdVi \JdVi 



1, 



< c 



Br{0)\V, 



\VvLl\^dV + 


([ \Vulfdv] 




\JBn(0)\P^ / 



1/2- 



when e is small enough. Now, combining (145|) and ( 146|) we obtain that 



1/2 

+ ( / \Vut\'dV^ 

\JBii{0)\Vi 



Wu'rdv 



L 



\ul\^ dV < 06^-"" 
Br{o)\v, UBR{o)\r, 

1/2 

\Jbr(0)\V^ 



Wu'rdV 



From P7|) . we find that 



Wu'rdV < Ce 



6-2a 



BR{0)\ri 



and then, from 



'BrWW, 

uniformly in £ ^ 1. Finally, fl57|) follows from ([35]) and 

Recall that R = e'^~°'. Rewriting ( I35ll37|) in terms of x gives 



/ 



B,^(^!)\Vf 



W^--N,YdV < Ce 



dVf 
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when e is sufficiently small. Furthermore 



9e / {\il,vf da = ge^ / {(u^ - w„ z/) + (w^, z/)}^ rfa 

JdVf JdVI 

/ {wi,iyfda + 2 (wj, z/)(uf - Wj, z/) (icr 

JdVf JdVf 



3-2a 



+ 



dVf 



By Holder's inequality, (|5Ti) . and the fact that w G C°°(f2), we have 



(wj, z/)(uf - Wj, z/) da 

\ 1/2 

< I / iwi.vfda] 

hvi J \JdVt 

1/2 



1/2 



[u^ - Wi, v) da 



< 



dVf 



1/2 



I uf — w,- 1 da 



dVf 



Since the last integral in (152!) is O (e^), we conclude that 

g, f (uf , da = ge'-'"' [ (w„ da + O (5^-2°) . 
JdVf Javf 

We now return to estimating [w^]. From flHT]) we have 

Vw^ = Vw + Vz^ , 

where 

Vz^ = 5^{0(e--|x-xf|)V(uf-w) 



X — X • 



X — X- 



(uf - w) 



Then, since the supports of </> (e |x — xf |) and |x — x^ |) are mu 
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nonintersecting for any i 7^ j G 1, . . . , A''^, using the definition of 0, we have 
/ \VztfdV<2j2 [ ^'{e--\^-x^])\Viul-w)fdV 
+ iE / [<P'{e-\^-^t\)T\^t-^\"dV (56) 



< 



|2 , 1 I _£ „ i2 



|V«-w)|^ + — K-w 



where C depends on only. Since w G C°°(n), the following estimates hold 



|uf (x) - w(x)|' < 2|uf (x) - w,|' + 2|w(x) - Wif 
<C[|u^(x)-w,|V|x-x^|2] , 
|V(u^(x)-w(x))|^<2|Vuf(x)|VC^ 



(57) 
(58) 



for every x G B^k (xf) \Vf, where C > is a constant that depends on w 
only. Therefore, by (H9l) and (l50ll we obtain 



\V{ut-w)fdV 



< 2 



+ C\Be^ (xf)| = o(e"^^'^^6-"'3''>) , 



(59) 



and 



< c 



Here I-B^k (xf)| is the volume of -B^k (x^). It follows that 



I uf — w,- 1 dV + e 



3k 



(60) 



Jn\UiVf 



(61) 



since 1<q;<2, l<K<a, and there are O {e ^) spheroidal particles. In 
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addition, by Holder's inequality, 

\\/w'fdV= [ \Vzl + VwfdV 

= I \VwfdV + o{l)= I dV + O {e^^''-^^) + o{l) 

= [ \VwfdV + oil), 
Jn 

(62) 

when 6 is small. This result extends to w G H^{Q) by a density argument. 

Next, consider the asymptotic behavior of the nonlinear term. Extending 
continuously w"^ to G H^{Q) and using the uniform boundedness of 
in H^{Q) (e.g. from ( |62|) and Poincare's inequality), we conclude that there 
is a subsequence such that ^ w weakly in H^{Q) and strongly in L^(f2) 
where 1 < p < 6. Since the Lebesgue measure of the set UjPf converges to 
zero when 5 — )■ and w G C°°{Q), we have that 

/ {l-\w'f)^dV^ [ {l-\wfYdV, (63) 

as £ — )■ 0. 

Finally, by f l54|) . we determine that 

ge'-^'^ V / (uf , da = ge^-^^ V / (w„ vf da + O (^^-2") 

= ge'-'"J2f i^u'^fda + oil) (64) 
= 9£'^Y] / {wi, i^f da + o{l) , 

A JdVi 



since a <2. Thus 

"^^^ dV 



£^[w^] = j J|Vw|'+ (1 - |w|') 

. J 9 



'11 



Wj, p)^ da + oil) 



(65) 
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when e is small. It remains to determine the asymptotic limit of the boundary 
term as e — )■ 0. The sum in this term can be rewritten as follows 



ge'' V / (w (x^) , uf da = ge^'Y^ f K (xf ) e^, uf da 
= 9£^y^ I {ek,v){ej,v)da 
= 9^'^y\ / {eu,Rlv){ej,Rlv)da 

, VJdV 

= [ (A^(x)w(x),w(x)) dV, 
Jn 



(66) 



X- ) Wk (xf) 



where e^, k = 1,2,3 is an orthonormal basis in 



the matrix-valued 



function A%{-x) = (/e'^ ^ . 5 (x — xf ) J^^, {e^, Rii'){ej, R^u) da, the function 

is a rotation matrix, such that Vi = RfV. Fur- 
ther, Wk = (w(x^),efc), where k = 1,2,3 and we assume summation over 
the repeated indices. Thus, from our assumptions on the geometry of the 
domain 



£e [w'] I |Vw|^ + (1 

for every w G C°°{^1). 



w 



W, WJ 



dV. 



(67) 
□ 



Theorem 4. Let a sequence of minimizers {u^} of be such that the se- 
quence {ug} of extensions o/{ue} to Q converges weakly in H^{Q) to some 
u e H^{n). Then 

liminf^^u^] > ^[u] , (68) 



e-s>0 



where E is defined by (E^. 



Proof. Suppose that there is {u^} C C^{fl) such that — )■ u strongly in 
(Q) and the extensions to Q of minimizers of Sg converge ^ u weakly 
in H^{Q). We construct u| = u^ + zf in the same way as in (13T|) . so that their 
extensions u| ^ converge weakly in H^{Q) along with Se [u|] — S [us] 
as £ — > 0. Let Q := Us — u| and denote its restriction to fl\ Uj V- by 
Then (^^ ^ (^5 := u — , weakly in H^{Q) and strongly in L^(f2) for p < 6 
as e — > 0. 

1. Asymptotics of(^. First, we show that 



lim lim sup VC, 



0. 



(69) 
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We begin by observing that the expression for ^^^[ug] can be rewritten so that 
£,[n,] = £,[nl] + [ \VC'S dV + 2 f (VCf , Vuf) dV 



+ / \C:\ dV-2 \C!\ dV-A (Cf,<) dV 

Jn\u^vf Jq\u,v^ Jn\\j,vi 



+ 2 [ \c!\'\nffdV + A[ {C!,4YdV (70) 
+ 4/ dV + A I \CWCul)dV 

+ 2Y,9s [ {C!,iy)'da + Y,9e [ (C',^) {niu) da. 

JdVf , JdVf 



Since is a minimizer of S^, we have that 
then 



/ \VC!\'dV<-2[ (VC^Vuf)d\/ 

+ 2! \Cl\'dV + A! {Ct.viDdV 

-4 / |ufr(C',u^) dV-A [ \C!\' dV 

-2Y,9e [ {C!,u)'da~Y.9e [ (C,^) da. 

JdVf JdVf 



(71) 



We need to estimate each term on the right hand side of (I7T]) . In the remain- 
der of the proof, C > denotes a constant independent of e and 5. 
(a). Beginning with the first term, we write 



iVulVCI)dV= / iVus,VCI)dV 



r (72) 

+ / {Vz^„VQ)dV. 

Jn\UiPf 
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We have 



/ 



Vus,VCs]dV- I (Vus,VCl)dV. (73) 



The second integral in f l75]) can be estimated with the help of Holder's and 
Minkowski's inequalities 



Vvis,VQ)dV< / \Vus\'dV 

1/2 



1/2 



1/2 



< 



1/2 



< c 



< c 



1/2 

\Vus\'dV] <C\ 

1/2 



1/2 



|Vu5 - Vu + Vul dy 



1/2" 



iVus - Vu\ dV 



1/2 



Vul'dV + / |Vu5- Vurdy 



1/2 



1/2 



when £ — )■ because u G if^(fi) and |UiPf | — )■ 0. Here C > is independeiit 
of 6. Consider now the first integral in ( 1731) . By the weak convergence of Q 
to U5 — u and Holder's inequality, we have that 

/ (Vus,VCl)dV^ [ {Vus,V{us-u))dV 
Jn ^ ' Jo, 

\ 1/2 / /• \ 1/2 

2\ / / ivv |2 



< C||u — U5||//i(Q) , 



when £ — )■ 0. 
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Now, for the second term in fl72l) . we find using (l^Tj) and the Holder's 
inequality 



(Vz|, VCD dv < 



1/2 



iVz^l 



n\u,rf 



1/2 



|VC 



^ 0. 



when e — )■ 0. 

(b). Consider the second term in f l7ip . Since Q converges weakly to 
u — in H^{Q) and strongly in Lp{Q) for p < 6 when e — )■ 0, we have that 



/ \C!fdV<[ 




^ dV ^ f \vL-VL5\^dV 




Jn 




5 



as e — i- 0. 

(c). Using Holder's inequality, we find that 



(C,u3 dV 



< 



\ut\^dV 



1/2 



dV 



1/2 



{n) , 



when £ — )■ 0. 

(d). Estimating in the same way as in (c), we obtain 



uff(Cf,uf) dV 



n\UiVi 



\ 3/4 

fdv) 



u — u^prfV") <C||u — u^l 



dV 



1/4 



by Sobolev embedding. 

(e). Estimating in the same way as in (c), we obtain 



ICl'(euf) dV 

1/4 



< 



\ 3/4 



dV^ 

\ 3/4 



W.t dV 



1/4 



(^J\usfdV^ (^J\u-usfdV^ <C||u-U5 



|3 
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by Sobolev embedding. 

(f). We use 021]) with A = 1 to obtain 



< C 



< c 



[ IVCI 


^dV+ [ 




icir 


dV 








[ IVCI 


Uv+ [ 


CI 


2 

dV 






Jn 









C \u- usf dV < C\\u- us\\m{n) , 
Jn 

by the strong convergence of u| and in L^^Q), 1 < p < 6, to and u, 
respectively. 

(g). Using the Holder's inequality, we get 

J29s [ (C', ^) (uf , z/) da 

^ E ( K \9s\ (c^ dv^ 

As in (f), applying ( 126!) with A = 1 we have that 

limsup / \ge\ {ul,uY dV < C / ju^l^rfV^, 

e-*-0 JdVf Jn 

then, with the help of (f) we obtain the estimate 
limsup y^de (Ce^^) i^t^) do- 

6^0 . JoVf 

<c(^j \vis\^dV^ (^j \u-us\^dVy < C||u-U5||^i(f,). 
Now, combining (a)-(f) and using (17T]) . leads to 

limsup f \VCt\^ dV <C\\M-us\\m 



(74) 
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when S is small. This, along with strong convergence of to u when (5 — )■ 0, 
proves the claim (1^ . 

2. Limiting behavior of Sslu^]. From flTOj) and fl7T|) we conclude that 



S,[u,] > S,[ui] + 2 / (VC', Vuf) dl^ 
-2 [ \Ct\'dV-Aj {ClviDdV 
+ 4 / \nff{C!,ul) dV + A [ |Cl'(C',u3 dV 
+ 2j29e [ {C!,iyfda + J29s [ (C',^) i^liy) da 

JdVf -■ JdVf 



(75) 



Thus 



lim inf [ue] > lim lim inf 6^ [uf] — lim lim sup / I VC^ I ^ dV 



and 

lim inf [u^] > £^ [u] , 

because lim^-^o [uf ] = ^^[u^] and £ is continuous with respect to the strong 
convergence of to u in H^{VL). □ 

From (EZ]) and (EHD it follows that ^[u] < ^[w] for every w € H^{9), 
hence u minimizes £ over H^{VL). 

4-2. Magnetic Energy 

Having established the asymptotics of the liquid crystalline component 
of the energy, we now turn our attention to magnetic interactions. Consider 
f ll2p for the prolate spheroidal particle V with semiaxes a > b and long axis 



oriented in the direction of z-axis. It is well known [16| that the solution to 
this problem in the exterior of V is given by 

(tanh~^(t) -t) 2, (76) 



(a2 - 62)3/2 



30 



in cylindrical coordinates {p,6,z), where t = ^ ^/^(a^ — 6^)^^^ and ^ is the 
largest root of 



p2 



Further, m is the density of the magnetic moment, so that m = ^^iSl^iDl^ qt^(\ 
z is a unit vector in the direction of z-axis. Assuming that a ^ (p^ + 2^) 
and expanding in a/ (p^ + z^)^^'^ , we find that 

+ O (2;(a/r) ) , (77) 



3r 



3 



where r = \J + z"^ = |x|. Note that the leading term in (I77|) is identical to 
that for a sphere of the same volume as V and centered at the origin 23 



The leading order term in the expansion of the magnetic filed H generated 
by the ferromagnetic particle V is given by 

H(x) = (^^x - zj + O ((a/r)5) , 

then 

|H(x)| = O {m{a/rf) , (78) 

when a/r <^1. 

Now consider the term corresponding to the magnetic interaction between 
the particles V- and P| for some i,j = 1 . . . , iVe. We have 



|mf| \m'\Vo\(V^)Vo\(Vf) 



(P / (79) 



then 



/ (m„ H,) dV = (iV2£6«+2/3i-3) ^ Q (^6a+2/3,-9^ ^ q ^ ^gQ) 

when e — !■ by our assumptions a and 

Finally, we consider the interaction between the external magnetic field 
and ferromagnetic particles. We have 



/ (m„h,)f/V = V / (mf,h,)rfr 

= [ (h,M") dV ^ [ (h,M) dV, 
Jn Jn 



(81) 
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by fll4p and 0151) . where M is the effective magnetic moment density. 

Combining the results for the hquid crystal and magnetic energies, we 
conclude that the minimizers of the family of functionals converge to a 
minimizer of the functional 

JTgfu] = / [|Vu|V (1 - |u|^)% (Au,u) -2(h,M)] rfl^, (82) 
concluding the proof of Theorem [TJ 

Remark 2. Suppose that the particles are distributed periodically in 
with their centers of mass positioned at the vortices of a cubic lattice with 
the side e. If we assume that there exists a continuous function i? : f2 — )■ 
Orth+ := {X e M3><3 . = J, detX = l} such that Rf = i? (xf ) for 

every i = 1 . . . and e > 0, then 



and 



M(x) = mi?(x)z, (83) 

A(x) = gR{x) {Jg^ V ® vda) -R^(x) 
= g /?(x)(Ai(z ® z) + A2 (/ - z ® z)))/?^(x) (84) 
= ^(Ai(M(x) ® M(x)) + A2 (/ - M(x) ® M(x)))) , 

where Ai and A2 are the two distinct eigenvalues of /^^ v®vda. The coupling 
terms in (!82|) then take the form 

^^hz2A (M, u)^ + ^lu^ - 2(h, M). (85) 

For a needle-like prolate spheroid with a high aspect ratio we have that 
Ai ^ A2 and the coefficient A := ^^^^^^"^^ in front of (M,u)^ has a sign 
opposite to that of g. Hence nematic molecules align perpendicular to M 
when A > and parallel to M when A < 0. Since the model in [3| assumes 
that |u| = 1, the middle term in ( 185|) can be neglected and the remaining 
interaction terms in (153]) coincide with those in ([T]) up to a difference in 
notation. 
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